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Abstract 

Wc consider the role of precision measurements of beta decays and light meson semi-leptonic 
decays in probing physics beyond the Standard Model in the LHC era. Wc describe all low-energy 
charged-current processes within and beyond the Standard Model using an effective field theory 
framework. We first discuss the theoretical hadronic input which in these precision tests plays 
a crucial role in setting the baseline for new physics searches. We then review the current and 
upcoming constraints on the various non-standard operators from the study of decay rates, spectra, 
and correlations in a broad array of light-quark systems. We finally discuss the interplay with LHC 
searches, both within models and in an effective theory approach. Our discussion illustrates the 
independent yet complementary nature of precision beta decay measurements as probes of new 
physics, showing them to be of continuing importance throughout the LHC era. 



1 Introduction 



Beta decays played a central role in determining the V — A structure of the weak interactions and in 
shaping what we now call the Standard Model (SM) [H [21 E]- We focus here on the set of semi-leptonic 
"charged current" (CC) processes that are mediated in the SM by tree-level W exchange, up to ra- 
diative corrections. In the SM the CC weak processes are characterized by two main features: (i) the 
hadronic and leptonic bilinear densities involved in the process have a dominant V — A component, 
with other types of couplings — V + A, S, P, T — arising at higher order in radiative corrections or in 
recoil momentum; (ii) the effective Fermi constants extracted in beta decays obey lepton universality 
as well as quark-lepton, or Cabibbo, universality, which is equivalent in the SM to the unitarity of the 
Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. Universality relations can only emerge once 
the process-dependent radiative corrections are removed. Currently precision beta-decay measurements 
involving neutrons, nuclei, and mesons are used to probe the existence of non-SM interactions which 
effectively induce violations of the universality relations and/or novel non-(\^ — A) structures or correc- 
tions to the dominant vector and axial-vector couplings. The low-energy charged-current-interaction 
Hamiltonian is sensitive to many classes of SM extensions. In this sense, beta decay measurements 
can be considered as "broad band" probes of physics beyond the Standard Model (BSM): while by 
themselves they do not allow us to reconstruct the ultraviolet dynamics, they provide, at 0.1%- level 
precision, powerful boundary conditions and diagnostics on virtually any TeV-scale SM extension. 

Considerable experimental progress is ongoing or expected in a few-year time scale on several fronts, 
using both cold and ultracold neutrons |il5l|6l|7l[Hll9l[ini[IIl[12l[l3l[Tll[l5l|16l[I3 trapped 
nuclei [211 ES] , and rare pion and kaon decays [231 EH ESI ES] . Some of the measurements plan to reach 
sensitivities between 10~^ and 10"'^; this makes such observables very interesting probes of new physics 
effects originating at the TeV scale, because such effects are expected to have size 0{{v/ A^smY), where 
V = (2v^Gp')~^/^ ^ 174 GeV and Absm denotes the mass scale where BSM physics appears. 

As in previous reviews [271 EHl EHl |30], the overall goal of this article is to discuss the discovery 
potential and discriminating power of planned precision beta-decay measurements with neutrons, nuclei, 
and mesons, in light of other existing precision electroweak tests and high-energy collider searches, 
such as at the Tevatron and the LHC. In order to achieve our goal, we work within an effective field 
theory (EFT) framework, in which the dynamical effects of new heavy BSM degrees of freedom are 
parameterized by local operators of dimension higher than four built with SM fields. All model-specific 
analyses of beta decays can be cast in the EFT language and the limits on the effective operators we 
derive can be readily converted into constraints on the parameters of any SM extension. In the absence of 
an emerging picture of new dynamics from collider searches, the EFT analysis is the first necessary step 
to establish the motivation and significance of this set of low-energy probes. Subsequently, we will also 
discuss well-motivated models such as the Left-Right Symmetric Model and supersymmetric extensions 
of the SM in order to show the discriminating power that combinations of beta decay measurements 
can have on explicit models. 

Probing short-distance BSM couplings through precision phenomenology of beta decays requires 
knowing the relevant hadronic and nuclear matrix elements to a precision comparable to the size of the 
new physics effects one could expect to appear. This means that one needs to know the hadronic matrix 
elements of the SM operators, that is, of the V and A currents, to the level of /Absm)^), i-e. of 10~^ 
or better. This is a necessary condition for beta decays to function as competitive probes: we are in 
search of a small BSM signal, and hence we need to know the SM "background" to a level comparable 

^In this review we consider the decays involving the hght quarks m, d, and s exchisively. 

^The EFT analysis can be applied to all low-energy probes of CC interactions. It is also valid for collider searches 
as long as the particles which mediate the new interactions are above particle-production threshold at the operating 
center-of-mass energy. In this case, a direct comparison of low-energy and collider constraints can be performed, as we 
discuss in Section ISl 
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to that of the signal for which we are looking. One also needs to know the matrix elements of the 
BSM operators, such as the S, T, P densities, because all the observables are sensitive to the product of 
the short-distance BSM coupling with the appropriate hadronic/nuclear matrix element. Consequently 
if a certain matrix element is suppressed, the sensitivity to the corresponding BSM coupling is also 
suppressed. Moreover, were such an anomalous suppression absent, the fractional uncertainty on the 
BSM matrix element still determines how well we can constrain that BSM coupling. For BSM operators, 
the precision required of the relevant hadronic matrix element is much less severe; an uncertainty at the 
C(10%) level is acceptable. Motivated by these considerations, we pay special attention to the hadronic 
and nuclear uncertainties which appear. 

This paper is organized as follows. In Section[2]we set up the theoretical framework for the analysis of 
all low-energy CC processes within and beyond the SM. In Section [3] we discuss the status of Cabibbo 



universality tests (Sec. 3.1) and lepton universality tests (Sec. 3.2) and explore the implications for 
BSM physics. In Section |4 we focus on differential decay distributions in beta decays and discuss the 
implications for non-(V^ — A) couplings. In Section [s] we explore the constraints on non-standard CC 
couplings that can be obtained from LHC data. In Section |6] we illustrate how the precision tests can be 
used to probe the parameter space of models such as the Left-Right Symmetric Model and the Minimal 
Supersymmetric Standard Model (MSSM), and we present our concluding remarks in Section [Tj 



2 Theoretical Framework 

2.1 Effective Lagrangian 

In this review we take the point of view that the Standard Model emerges as the low-energy limit of a 
more fundamental theory characterized by the scale A at which new particles appear. Consequently, at 
energies scales below A, namely, A > E > Mz,Wi the new degrees of freedom are no longer present; they 
have been "integrated out," yielding an effective Lagrangian comprised of the SM Lagrangian augmented 
by a string of d > 4 operators constructed with the low-energy SM fields, suppressed by A'^~^ [HT] , that 
respect the SU(3)cxSU(2)iXU(l)y gauge symmetry of the SM. Flavor physics observables constrain 
the appearance of non-SM invariant operators to energies far beyond the weak scale [321 |33l |3l]. The 
building blocks of the gauge-invariant local operators are: the gauge fields G^, W^, B^, corresponding 
to SU(3)cX SU(2)iX U(l)y, respectively, the six fermionic gauge multiplets, including a singlet right- 
handed neutrino state, 

\ „,i „,i ji ji li / 



the Higgs doublet (p 



and the covariant derivative 



^0 ) , (2.2) 



A 



A 



D, = Id, - igs^G^, - ig^W; - ic^YB, . (2.3) 

In the above expression, / is the identity matrix, A^ are the SU(3) Gell-Mann matrices; a" are the SU(2) 
Pauh matrices; gg, g, and g' are the gauge couplings; and Y is the hypercharge of a given multiplet. The 
introduction of three light, right-handed neutrinos to accommodate the existence of neutrino masses 
illustrates explicitly that new light degrees of freedom can be included in the EFT if we assign their 
SU(3)cxSU(2)LxU(l)y quantum numbers. The impact of other light, new physics is left as an exercise 
for future work. 
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The leading operators which modify CC interactions are of dimension six, though it is worth noting 
that new physics can also modify the radiative corrections and hence the couplings with which the 
SM operators appear. The minimal set of operators contributing to low-energy semi-leptonic charged 
current processes can be divided into two groups: operators built out of SM fields, noting the left 
column below, in which we follow the notation of Refs. [351 EH], and operators involving the singlet 
R-handed neutrino field v [37], which are displayed in the right column below. Furthermore, within 
each group the operators can be divided into two classes — four-fermion contact interactions and vertex 
corrections. The vertex correction operators are written in SU(2)-invariant form and therefore involve 
the Higgs doublet: after electroweak symmetry breaking (EWSB) they include terms involving a W (or 
Z) boson, a fermion, and an anti-fermion. Here is the list: 
Four-fermion operators: 



off = {l^^an){q^^a\) 
Oiedq = {le){dq) +h.c. 
Og. = {lae)e^\qbu) + h.c. 
0Ll = (^aa'^'^e)e»^fea^.«)+h.c. 

Vertex corrections : 

O^ud = i{'^^(^D^ip){uYd) + h.c. 
Og=(y.t,B;y,)(g7Xg) 



Oevud = {e-i^v){u-i^d) + h.c. 
Oquv = {lv){uq) +h.c. 

'luqd 



O^L = {layY\qbd)+\i.c. 



r,(3) 
^luqd 



{qhcr^ud) + h.c. 



O' 



ipue 



i{(p^eD^if){v'^^e) + h.c. 



(2.4a) 
(2.4b) 
(2.4c) 
(2.4d) 

(2.5a) 
(2.5b) 
(2.5c) 



Denoting by Aj the effective dimensionful coupling associated with the operator Oj, we can write 
the effective Lagrangian as 



^(eff) ^ r 



SM 



1)2 



«i Oi 



with d,- 



V 



(2.6) 



where in the last step we have set the correct dimensions using the Higgs vacuum expectation value 
(VEV) V = {(f^) = {2\/2Gf)~^^'^ and defined the dimensionless new-physics couplings dj, which in 
general are matrices in both the quark and lepton flavor spaces. In this framework one can derive the 
low-energy effective Lagrangian at 0{1 GeV) for semi-leptonic transitions. It receives contributions 
from both IV-exchange diagrams, with modified I^-fermion couplings, and the four-fermion operators. 
After including the electroweak radiative corrections to the SM operator [3H], the matching procedure 
leads to a low-energy quark level effective Lagrangian involving ten dimension-six operators: 



CC 



Gp^Vud 



V2 



(^1 + 6,3^ 67^(1 - 75)z/e ■ M7''(l - 75)C? 



(2.7) 



+ 67^(1 - 75)z/£ ■ M7^(l - 75)c?-h eL e7^(l 75)//^ ■ M7''(1 - 75)c? 

+ e/j e-ff,{l - -fr,)ue ■ u-f^{l + -f5)d + in 67^,(1 75)^'^ ■ m7''(1 + 75)*^ 

+ es e{l - -f5)ui ■ ud + is e{l + ■ ud 

- tp e(l -75)z/£-M75(i - Ip e(l-F75)z/^-M75rf 

+ ct ea^u{l - 75)z/^ ■ ua^''{l - 75)6/ + et ea^y{l + 75)//^ ■ ua^''{l + 75)^ + h.c. . 

In the above equation cf^ /\/2 = g'^/{8M^) is the tree-level SM Fermi constant, and Sjs encodes the 
effect of SM electroweak radiative corrections to semi-leptonic transitions, noting that the Fermi theory 
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QED contributions have been subtracted [391 HOI EHl SH SSI US]. The couphng Gp^ can be expressed 
in terms of the Fermi constant = 1.166371(6) x 10~^GeV~^ precisely measured in muon decay 
In order to do so, one has to consider the low-energy effective Lagrangian describing muon decay 

^t^^ep^u^ = -4:GP {l + 6f, + e^) e£7^i/eL ■ '^(.lYi^'l + h-c , (2.8) 

where = G'^''(l + 6^ + e^). Here 5^ represents the SM electroweak radiative corrections [JB] to 
purely leptonic transitions, noting that the Fermi theory QED contributions have been subtracted, and 
encodes possible new physics contributions, so that Gp = G^{1 — — e^) 



The BSM effective couplings in Eq. (2.7) are denoted by and e/3, using the self-explanatory 



notation a,(3 = L, R, S, P,T. These couplings can be expressed in terms of the weak-scale couplings 



(Xj im mi EZ]. In the effective Lagrangian of Eq. (2.7), e,u, and d denote the electron, up-, and 
down-quark mass eigenfields, whereas ui represents the neutrino flavor fields. In general we can have 
i 7^ e — in what follows, we suppress lepton flavor indices. Finally, identical CC effective operators 
appear for other quark flavors. For example, the operators obtained by replacing the d quark with the 
strange quark s describe [AS"! = 1 semileptonic processes. 



Next, we discuss some noteworthy points in regards to the effective Lagrangian of Eq. (2.7): 



The effective couplings denoted by involve L-handed neutrinos, whereas e/3 involve R-handed 
neutrinos. Therefore, the appear in decay rates and distributions either quadratically or linearly, 
but the latter appears multiplied by the small factor rriy/Ey, as it is realized through interference 
of the SM and BSM couplings. In constrast, the couplings contribute linearly to the decay 
rates without niy/Ey suppression. As a consequence, the bounds on the e's are much stronger 
than the bounds on the e's. 

There are twelve SU(2)LxU(l)-invariant operators that contribute to beta decays, though there 
are only ten quark-level U(l)£;Af -invariant operators. This is because the correction to the SM 



operator encodes contributions from three weak-scale operators of Eqs. (2.4) and (2.5), namely, 
the contact operator Oi^ and the quark and lepton vertex corrections, Olpq and O^/. All other 
low-energy operators are in one-to-one correspondence with the TeV scale SU(2)2,xU(l)-invariant 
operators. It is interesting to note that SU(2) gauge invariance implies that the same couplings 
mediate not only charged-current processes but also "neutral current" processes such as ee -H- 
MM, dd. 

While the physical amplitudes are renormalization scale and scheme independent, the individual 
effective couplings es,p,T {^s,p,t) and the corresponding hadronic matrix elements display a strong 
scale dependence in quantum chromodynamics (QCD) already at one-loop order (see Ref. |18] 
and references therein). Throughout the paper, we quote estimates and bounds for the (e^) at 
the renormalization scale /z = 2 GeV in the MS scheme, unless otherwise specified. 



The Lagrangian of Eq. (2.7) mediates all leading, low-energy charged-current weak processes in- 
volving up and down quarks. In some charged-current processes involving first-generation quarks the 
theoretical and experimental precision has reached, or will soon reach, a level that allows stringent 
bounds on new-physics effective couplings. To set the stage for this discussion, we now provide an 



overview of how the various BSM couplings of Eq. (2.7) can be probed experimentally — we explore 
these points in detail in the following sections. For context, we note that detailed expressions of the 
noTi-iy — A) contributions to neutron and nuclear beta decay correlation coefficients can be found in the 
papers by Jackson, Treiman, and Wyld IHI] , where one can re-express the Lee- Yang couplings 



they employ in terms of the and using the expressions given in Eqs. (2.17) below. 



^Our notation in Eqs. (2.7) and (2.8) corresponds to that of Ref. [30] if we replace 5p A-fj^ and 5^ Ar 
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The combinations [e^ ± e^j) affect the overall normalization of the effective Fermi constant in 
processes mediated by the vector and axial- vector current, respectively. As discussed below, the 
hadronic matrix elements of the vector current are known very precisely up to corrections due to 
QCD flavor symmetry breaking, that is, quark mass differences, whereas the axial-vector matrix 
elements require non-perturbative calculations. Therefore, while the difference (e^ — en) remains 
relatively unconstrained, the sum {ei + Gr) is strongly constrained by quark-lepton universality 
tests, which are tantamount to CKM unitarity tests. These tests involve a precise determination 
of Vud and Vus from processes mediated by the vector current, such as 0^ — )■ 0"*" nuclear decays 
and K — )■ Triu. An extensive analysis of the constraints on {ei + ej?) from universality tests 
and precision electroweak observables at the Z-pole was performed in Ref. [45J. In this context 
it was shown that constraints from low-energy are at the same level or stronger — depending on 
the operator — than those from Z-pole observables and e~^e~ — )■ qq cross-section measurements at 
LEP. 

The right-handed coupling en affects the relative normalization of the axial and vector currents. In 
all beta decays en can be absorbed in a redefinition of the axial coupling, and, up to calculable ra- 
diative corrections |l2l [52l [531 [511 ESI [56] , experiments determine the combination {l — 2efi)gA/gv, 
where gy and gA are the vector and axial form factors at zero momentum transfer, to be precisely 
defined below. Disentangling requires precision measurements of (1 — 2eji)gA/gv and precision 
calculations of gA/ gv in lattice QCD, which, unfortunately, are not yet at the required sub-percent 
level. 

The effective pseudoscalar coupling ep contributes to the leptonic decays of the pion. It is strongly 
constrained by the hehcity-suppressed ratio R^, = T{7i — )■ ei'['y])/T{7T — )■ /ii/[7]). Moreover, as 
discussed in Refs. ESI EH], the low-energy coupling ep receives contributions proportional to 
es,T through electroweak radiative corrections. 

Both the scalar and tensor couplings and contribute at linear order to the Fierz interference 
term b in the beta decays of neutrons and nuclei, as well as to the neutrino-asymmetry correlation 
coefficient B in polarized neutron and nuclear decays. The empirical determination of the beta- 
asymmetry correlation coefficient A and the electron-neutrino correlation a in neutron and nuclear 
beta decays, as well as positron polarization measurements therein, entrain sensitivity to the Fierz 
interference term as well. Thus bounds on es and can also be obtained from these observables. 
Moreover, the quadratic dependence on these couplings is useful in limiting their imaginary parts 
as well. Finally, the tensor coupling can also be constrained through Dalitz-plot studies of the 
radiative pion decay vr — )■ ez/7. 

Neglecting neutrino masses, all the couplings contribute to decay rates as per oc |e^P, so that 
it is more challenging to set limits on their appearance at low energies. 



• All of the operators of Eq. (2.7) can produce collider signatures. Before the advent of the LHC, 
collider bounds on the chirality-flipping scalar and tensor couplings es^p^r and is^p^r were very 
weak, because their interference with the SM amplitude appears with factors of rrif/Ef, where rrif 
is a light fermion mass with / G {e, u, d}, which at collider energies strongly suppresses the whole 
effect. At the LHC, however, the contributions which appear as |e/3p or |e/3p can be boosted by 
a factor involving the energy in the numerator, noting that we replace (w/Absm)^ (-^/Absm)^, 
thus increasing the sensitivity to these couplings. We will discuss these bounds and show that with 
higher center-of-mass energy and integrated luminosity they become competitive with low-energy 
searches for es,T or stronger than low-energy bounds for ep^s,T- This analysis, of course, makes 
sense only for Absm ~ few TeV. 

The above considerations and more are summarized in Table [U 
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Table 1: Summary of the most sensitive direct low-energy probes of non-standard cliarged- 
current couplings. Left column: combination of couplings. Right column: probe. The 



effective couplings are defined in Eq. (2.7). The decay parameters a,b,B,A are defined in 



Eq. (4.41). If the new interactions originate at mass scales above the TeV, the LHC provides 



constraints on all non-standard couplings through the process pp ^ e + u + X. 



Non-standard coupling 


Probe 


cl + ep 


j^'^ X /■ TV T 'J 'J 

CKM unitarity 


— e_R, cp, cp 






0, B [a, A\ 




b, B [a, A], 7r — )■ ez/7 




a, 6, 5, A 



2.2 Hadronic and nuclear matrix elements 

Hadronic and nuclear transition amplitudes always involve products of short-distance couplings, evolved 
to the appropriate matching scale, and hadronic matrix elements. Thus in order to extract information 
on the former, we need to know the latter. Specifically, we need to match the quark-level effective theory 



of Eq. (2.7) to a low-energy effective theory written in terms of meson and baryon degrees of freedom. In 
QCD, this effective theory is Chiral Perturbation Theory (ChPT) [201 EH 112] • In the baryon sector, the 
low-energy structure of the theory is more complicated, and heavy baryon chiral perturbation theory 
is employed [63], where we refer the reader to Ref. [6l] for a review. Different systematic approaches 
to remedy its limitations have been developed, improving the theory's convergence, notably the "small 
scale expansion" of Refs. [651 ESI Ell |68], as well as Ref. [69]. As we have discussed, the precision with 
which we know the matrix elements of the SM operators limits our ability to constrain new physics. 
If we wish to probe scales such that ("w/Absm)^ ~ 10~^, we need to know the SM matrix elements 
with commensurate precision. This requires including all of the electromagnetic, isospin-breaking, and 



recoil-order effects in the calculation. Since the operators appearing in Eq. (2.7) have the factorized 
structure Jquark X </iepton, we need not present the ChPT framework but rather can describe the purely 
hadronic effects in terms of meson and nucleon matrix elements of quark bilinears. Nevertheless, the 
full ChPT machinery should ultimately be employed to compute long-distance radiative corrections. In 
the case of neutron decay, this has been done in Ref. [70], finding results consistent with non-ChPT 
based calculations [121 IH21 ESI EH EHl EHl [71]. In this review, we will not further discuss long-distance 
radiative corrections to neutron decay and refer the reader to Refs. [TOj and j42j for recent detailed 
accounts. 
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2.2.1 Meson matrix elements 



Leptonic (M — )■ lu) and semi-leptonic (Mi — M2I1') decays of pseudoscalar mesons provide strong 
constraints on the CC BSM couplings. The relevant one-meson matrix elements are parameterized in 
terms of the pion and kaon decay constants Ft^^k as follows (in our normalization F^^ ~ 92 MeV): 

(0|M7M75rfk-(p)) = -iV2F^p^ (2.9) 

{0\uj^-f5s\K-{p)) = -iV2FkP^ (2.10) 

2 

(0|M75rfk'(p)) = I "^^ (2.11) 

2 

{Q\u-i,s\K-{p)) = I V2Fk. (2.12) 

ruu + rus 

The pseudoscalar matrix elements follow from the axial-vector ones by using the operator relation 
Qil^l^lj = iiiT^i + fTij) Qjl^Qi- Matrix elements of the other bilinears vanish by parity. 

The two-meson matrix elements of vector and scalar densities can be parameterized in terms of 
two form factors. Specializing to the — )■ 7r~ transitions, we have, noting the momentum transfer 
q = p — k: 

{n-{k)\-s^,u\K\p)) = {p + k)^U{q') + {p-k),f.{q') (2.13) 

2 2 

'Yji ^ nnrt ^ 

1,7:- {k)\MK\v)) = -fo{q'), (2.14) 

nis - rriu 

where foiq"^) = f+{q^) + {(1^ /(j^k ~ ^n))f~il'^) the operator relation dfj^qi'-y^qj = i{mi — rrij) qjqi 
has been used to relate the vector and scalar matrix elements. In the SU(3) / limit the light quark 
masses obey = rrid = rris, so that /+(0) = 1 — the Ademollo-Gatto theorem [721 [73] ensures that the 
corrections to the flavor symmetry limit start at second order: /+(0) = 1 + 0{{ms — m^)^)]^ Finally, 
for completeness, we report the tensor matrix element, which involves a new dynamical form factor 

{n-{k)\7sa,MK'{p)) = I B^iq') . (2.15) 

niK 

The decay constants and form-factors can be calculated in lattice QCD (LQCD) and we will review 
the relevant results as needed. 



2.2.2 Nucleon matrix elements 

At the one-nucleon level, we require the matrix elements between the neutron and proton of all possible 
quark bilinears of dimension three. These can be parameterized in terms of Lorentz-invariant form 

''Note, however, that in the case of charged Kf^ decay the existence of 7r° — 77 , 77' mixing imphes that corrections to 
/+(0) = 1 occur at first order in {irid — mu)/{ms — m), with m = (md + m„)/2. 
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factors as follows [75]: 



{p{Pp)\ ujf,d\n{pn)) 

{p{Pp)\ u-rf,-f5d\n{pn)) 

{p{Pp)\ ud\n{pn)) 
{p{Pp)\ u-f5d\n{pn)) 

{p{Pp)\ua^^d\n{pn)) 



UpiPp) 
Up{Pp) 



UniPn 



9p{q^ 



2Mn '^^"^ ^ 2Mn 



(2.16a) 
(2.16b) 



= 9s{q^) Up{pp)Un{Pn) (2.16c) 

= 9p{q'^) Up{Pp)'j5Un{Pn) (2.16d) 

= Up{Pp) 9T{q^)(^t,u + 9?{q'^){qtiiu - quip,) 

+ 9?\q^) {q^Pv - quP^) + 9T\q^) [i4iv - iviii^j Un{Pn) , (2.i6e) 



where Up^n are the proton and neutron spinors, P = p^ + Pp, q = Pn — Pp is the momentum transfer, 
and Mjv = (M„ + Mp)/2 denotes a common nucleon massF] Note that the above spinor contractions 
are C(l), except for Up^y^Un, which is O^qjMj^). 

In order to make contact with the standard references on neutron and nuclear beta-decay phe- 



nomenology [511 HSl EHl EH], we note that upon neglecting recoil order terms Eq. (2.16) can be viewed 
as the matching conditions from our quark- level effective theory Eq. ( 2.7[ ) to the nucleon- level effec- 
tive theory originally written down by Lee and Yang [3T]. The Lee- Yang effective couplings Cj, C[ 
[i e {V, A, S, T}) can be expressed in terms of our parameters as follows P3| 



Cv 
C'v 
Ca 
C'a 
Cs 

C's 
Cp 

C'P 



G 



(0) 



Vv A Ci 



72 

gv {l + 5ii + eL + eR + eL + Ir) 
gv {l + Sfj + ei + eR-iL- Ir) 

-gA (1 + 5/3 + cl - Er - Cl + Er) 

-gA (1 + 5/3 + - + - Ir) 

9s (es + £5) 

9s {^s - £5) 

gp (ep - ep) 

gp (ep + ep) 

4 gx {er + e^) 

igr {er - It) ■ 



(2.17a) 

(2.17b) 
(2.17c) 
(2.17d) 
(2.17e) 
(2.17f) 
(2.17g) 
(2.17h) 
(2.17i) 
(2.17j) 
(2.17k) 



Using these relations and the results of Ref. |19] one can easily work out the dependence of beta decay 
observables on the short-distance parameters e^ and e^. 

Our goal is to identify TeV-scale induced new physics contaminations of typical size ea ~ (f /Absm)^ ~ 
0(10^'^) to the decay amplitude, so that they are comparable in size to the recoil corrections of 
0{q/Miy) ~ 10^^ and the radiative corrections of O^a/n). Thus, it is useful to organize the discussion 
in terms of a simultaneous expansion in new physics contributions, recoil, and radiative corrections 



^In the case of vector and axial bilinears, the Gordon decomposition can be used to trade the induced tensor term 
proportional to a^j,^q'^ for an independent scalar term proportional to P^. Here we choose to follow the parameterization 
of Ref. 

^ Various metrics and conventions appear in the literature. Lee and Yang [51j employ the "ict" metric, which in this 
case maps to the metric we employ if we let 75 — >■ —75 in their effective theory, noting 75 — 17*^7^7^7'^. Refs. |76| and 
pi] employ the metric we do but flip the sign of the 75 terms. 
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keeping terms through first order only. Higher-order terms may prove neghgible in hght of anticipated 
experimental sensitivities, but we indicate the role of certain, more significant ones. Employing this 
simultaneous expansion in e^, q/Mj^, and a/vr, we now discuss the contributions from the quark-bilinear 
operators: 

• Vector current: The form factor 5'y(0) contributes at 0{1) to the amplitude, whereas ^t(v)(0) 
and ^5(0) contribute at 0{q/Mi\f). Also, in the SU(2)j, or isospin, limit, the weak magnetism form 
factor ^T{y)(0) is related to the difference of the empirical proton and neutron magnetic moments, 
which are well-known, and the induced-scalar form factor gsiq"^), reflective of the presence of a 
second-class current, vanishes Corrections to the isospin limit are of C((M„ — Mp)/M]^) ~ 
q/Mjy. Since gs multiplies one power of q^/Mjy, its contribution to the decay amplitude is effec- 
tively of second order in the recoil expansion. 

• Axial current: The form factor 5'a(0) contributes at 0{1), whereas ^t(A)(0) and ^p(O) contribute 
at 0{q/MN). The induced-tensor form factor gT{A){q^) vanishes in the isospin limit [75], so that 
its contribution to the decay amplitude is of second order in qjM^. Similarly, the contribution 
associated with the induced-pseudoscalar form factor is quadratic in our counting, because the 
pseudoscalar bihnear is itself of order qjM^, and is accompanied by an explicit q/M^ suppression 

it can be studied in muon capture, note Ref. [771 EH] for reviews. 

• Pseudoscalar bilinear: The pseudoscalar bilinear Up'j^Un is itself of order q/M^q. Since it 
necessarily multiplies a BSM effective coupling ep because there is no pseudoscalar coupling in 
the SM, this term is also of second order in our expansion. 

• Scalar and tensor bilinears: These bilinears enter into the analysis multiplied by new-physics 
effective couplings ts,T- Computing the corresponding matrix elements to zeroth order in the recoil 
expansion suffices to identify (73(0) and gri'^)- Note that gT''^'^\q^) appear only in 0{q/Mi\f), and 
g^\q^) vanishes in the isospin limit [75j . 

In summary, to the order we are working, the amplitudes depend only on gi = gi{0) with i G 
{V, A, S,T} and ^T(y)(0). Up to second-order corrections in isospin breaking, one has gy = ^ [7 ^ 17 ^ 
I5U] . We define the ratio of the axial to vector form factors as A = gA/gvi where A > under our 
conventions. As we have noted, the neutron-decay amplitude in the presence of non-standard right- 
handed interactions is actually a function of A = A(l — 2eR). The parameter A is extracted very precisely 
from beta-asymmetry measurements in polarized neutron decay, leading to A = 1.2701(25) [81]; this 
number is essentially gA- There are no direct experimental handles on gs,p,T- 

A first principles calculation of gA,s,p,T, however, is possible with LQCD. The status of LQCD 
calculations of these charges is critically reviewed in [17], and in which the first estimate of gs from 
LQCD is provided. Different calculations give results in the range 1.12 < gA < 1.26; the errors are 
much larger than the experimental uncertainty. In constrast, the new estimates for the scalar and tensor 
charge in the MS scheme and at /i = 2 GeV are gs = 0.8 ± 0.4 and gT = 1-05 ± 0.35. Besides statistical 
uncertainties, which are particularly large for gs, the dominant LQCD systematic effects in gs^r arise 
from extrapolation in the quark mass to the physical point and from the renormalization constants, 
noting perturbative calculations were used to arrive at the reported results, and the non-perturbative 
calculation is in progress. More recently Ref. [82] provided improved results for the scalar and tensor 
charges, gs = 1-08 ± 0.28 and gx = 1.038 ± 0.011, with the uncertainty associated with statistics 
and chiral extrapolation. These results, however, do not include an estimate of the systematic error 

'^This term, however, is enhanced. Using the partially conserved axial current one can show that the form factor is 
of order M^r/rTig ^ 100, making a 0(10"^) contribution to the amplitude. The effect of gp on beta-decay rates has been 
worked out in Ref. j76j , and it should be included when the experiments reach that level of precision. 
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associated with finite volume and finite lattice spacing extrapolations. Therefore, in what follows we 
use the results of Ref. [?7] as the basehne lattice results. 



2.2.3 Nuclear matrix elements 



In moving from neutron decay to a general nuclear beta decay, a number of important differences 
appear. For one thing, numerous spin sequences are possible. Also the daughter state may well itself be 
unstable under electromagnetic or strong interactions. Finally, the Q-value, defined as the maximum 
electron (positron) kinetic energy, is generally much larger than the 0.8 MeV found in neutron beta 
decay and in some cases can be as large as 10-15 MeV. The electron/positron energy dependence of the 
decay observables are then much more appreciable. 

The spins are easily dealt with by the use of Clebsch-Gordan coefficients Cjf j^j'*^ together with 
the Wigner-Eckart theorem to define reduced matrix elements. Thus the general form for a vector and 
axial- vector nuclear matrix element between parent and daughter states having spins J, M and J', M' 
and masses Mi and M2, respectively, is \76y\ 



< /3\Vf,\a > 



P ■ £ 



(1 



2Ma 



2Ma 



.h{q 



Sjj'hiM' + ii^YJ-Cf,^^'f {q X £)fc 



2Ma 



. f^M'k;M 



-^ll;2 ^nqn' + ^^^^^ 



2Ma 



qn' + 



:P 



Air 




AMa 



(2Ma)2 



h{q^)q^P'^q 



{2Ma) 



-Mq') 



. ^M'k;M^nn';kr) /^^vn'/^x ^ 



{2MaY 



(2.18) 



where i'^ denotes the leptonic current and Ma = {Mi + M2)/2. Here each term corresponds to one in 
the analogous neutron transition via 

^ Qv, c qaV^ 

b^gT{v)VS, d-^gT(A)V^ (2.19) 
e^gs, h-^gpy/3. 

In addition, there exist terms f,g,j2,j3 which have no J = ^ — t- J' = | analog since they involve 
AJ = 2,3. 

For each form factor there exist known one-body operator, or impulse approximation, predictions. 



®Here we discuss only allowed decays, for which AJ = 0, ±1 with no change in nuclear parity. 
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Defining the nuclear mass difference A — Mi — M2, we liave 



b{q^) ~ A[gTiv){q^)MGT + 9v{q'')ML] 

c{q') ^ {l+^r'9A{q')[MGT+l{q'-A')M.r. + -^Miy{2A' + q') 



d{q') ^ {l + ^)-W)[-MGT-l{q'-A')M.r^ 

+ AM„L + MAM^rp] ± AgT(A) {q^)MGT 

e{q') ^ (1 + ^)~'9v{q')[MF + ^(g' - A')Mr^ - ^Mr-p] ± A~gs{q')MF 

f{q') ~ gv{q')2MAM{r,p} 

9{q') ^ -9v{q')^MQ 

h{q') ^ -{l+^)-'x[9A{q')-^Miy + Mq')A'MGT] 



ji{q') ^ -^9Aiq')Miy, i = 2,3, (2.20) 



wliere tlie upper (lower) sign refers to electron (positron) emission and the A^'s represent reduced, 
nonrelativistic matrix elements, namely. 



Mf = </3|lE- 

i 

Mgt = </3|lE' 

i 
i 
i 

Mr.p = ^<P\\J2^Hri-Pi+Pi-ri)\\c^> 

i 

MaL = < (^W^T^ai X {ri X pi)\\a > 

i 

MKy = ^<^\\Y.^rrlctKXC7,r.Y^'{h)\W>, X = 1,2,3 

i 

^^rp = ;t^<^I|5]^^[{^^-^"'-^"'^} + {^^-^"'-^"^}]II«>' (2-21) 

ZlVl AT 
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as well as 



Mq = <(3\\J2^MY,\h)\\a> 

i 

^ i 

Ml = < J] rf(rl X > . (2.22) 

i 

One expects meson exchange, or two-body, corrections to these predictions at the order of 5 to 10% or 
so 

If we neglect recoil effects, which is generally a good approximation because they enter the nuclear 
matrix elements at 0{A/Mn), there exist only the leading Fermi and Gamow- Teller form factors a = 
gv-Mp and c = qa-Mct- Most experiments are analyzed in terms of only these quantities. In the SM 
the vector weak charge is also the isospin raising operator, so that the Fermi form factor a(0) must 
vanish unless the parent and daughter states are isotopic analogs. Moreover, if the parent-daughter 
analog states have spin-parity 0"*", then the Gamow- Teller form factor vanishes — c(g^) = — and the 
Fermi matrix element becomes a simple numerical factor, namely, for unit isotopic spin a(0) = a/2, 
which is an exact prediction up to isospin breaking effects. The calculation of the latter effects has been 
an area of active investigation (HH [851 ESI EH EHl ESI [90l [91]. In simple terms, one needs to take into 
account that the "last" proton in the parent positron emitter is less strongly bound than the "last" 
neutron in the daughter state. An extensive analysis of data by Hardy and Towner, taking into account 
isospin breaking, has yielded the very precise value Vud = 0.97425(22) |92] . 

No symmetry principle determines the size of Gamow- Teller matrix elements, and these must either 
be determined empirically from lifetime or correlation coefficient measurements or calculated from 
nuclear wavef unctions, though the latter can only be done precisely if meson exchange effects are 
included. In the case of BSM matrix elements, wavefunction calculations of the requisite matrix elements 
are required. 



3 Cabibbo and lepton universality 

Beta decays provide stringent constraints on non-standard couplings through two classes of observables. 
(i) Total decay rates, after inclusion of radiative corrections, provide information on the strength of 
weak interactions, thus enabling precision tests of Cabibbo and lepton universality, (ii) Differential 
decay distributions, including spectra and correlations, are sensitive to the Lorentz structure of the 
underlying weak interaction, thus enabling searches for small non-(\^ — A) components. In this section 
we focus on universality tests. 



3.1 Cabibbo Universality 



In the SM the effective Fermi constant controlling semi-leptonic transitions Ui ^ dj for i,j G 1, 2, 3 
is related to by 

Gi3 = G^V,j (l + 5/3-(5^), (3.23) 



where G 



/9 



(0), 



VijGp 



1 + Sjs), the unitary CKM matrix Vij parameterizes quark mixing, and and 



Sfj, encode electroweak radiative corrections (see Eqs. (2.7) and (2.8)). We can thus test Cabibbo (or 



quark- lepton) universality by testing whether \Vud\ + \Vus\ + \Vub\ = 1- Beta decay rates permit 
access to the CKM matrix elements Vud and Vus- Since both the SM prediction and the experimental 
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measurements have reached the sub-percent level, these observables provide strong constraints on new 
physics, through the parameter Ackm defined as 



AcKM = \Vudr + \Vusr + iVubl' - 1 • (3.24) 

Here Vij are the CKM elements determined phenomenologically from semileptonic decays assuming 
only SM dynamics. The ratio Vij/Vij is parameterized in terms of the BSM couplings, as exemplified 
by Eqs. (3.25) and (3.29). In the unitarity sum, \Vub\ = 3.51(ii4) x 10~^ [HI] plays no role; Vud and Vus 



are both important and can be determined with high precision in a number of channels. The degree 
of needed theoretical input varies, depending on the component of the weak current which contributes 
to the hadronic matrix element. Roughly speaking, one can group the channels leading to Vud,us into 
three classes: 

• Semileptonic decays in which only the vector component of the weak current contributes: These 
are theoretically favorable in the SM because the matrix elements of the vector current at zero 
momentum transfer are known in the SU(2) j (SU(3) f) limit of equal light quark masses. Moreover, 
corrections to the symmetry limit are quadratic in mg^d — [ZSl [Z2] • Super-allowed nuclear beta 
decays (0"*^ — )■ 0"^), pion beta decay (tt^ — )• n'^e'^Ue), and K — )■ tt^z/ decays belong to this class. The 
determination of Vud,us from these modes requires theoretical input on radiative corrections [HI 
Sa Sal [931 El [95] and hadronic matrix elements via analytic methods |96l [971 |98l [HHl HOOl [Ml |I02] , 
or lattice QCD methods [1031 [M HOSl [T06l[T07] . 

• Semileptonic transitions in which both the vector and axial component of the weak current con- 
tribute: Neutron decay [n — )■ peu) and hyperon decays (A — >■ pep, ....), as well as nuclear mirror 
transitions, belong to this class. In this case the matrix elements of the axial current must be 
determined experimentally |108[ 1109] . Inclusive r lepton decays r — t- hur also belong to this 
class as both V and A currents contribute, though the relevant matrix elements can be calculated 
theoretically via the Operator Product Expansion 1111] . 

• Leptonic transitions in which only the axial component of the weak current contributes: In this 
class one finds meson decays such as 7i{K) — )■ fiu but also exclusive r decays such as r — )• Urir^K). 
Experimentally one can determine the products VudF^ and VusFk- With the advent of precision 
calculations of Fx /Ft, in lattice QCD |1121 11131 1114[ 11151 1116j . this class of decays provides a 
useful constraint on the ratio Vus /Vud jllTj 1118] . 



Currently, the determination of Vud is dominated by 0^ — )• 0+ super-allowed nuclear beta decays [96l 
leading to Vud = 0.97425(22), while the best determination of Vus arises from K£2,K£3 decays, 
leading to the best fit Vus = 0.2256(9) [miinH]. These determinations lead to Ackm = (1 ±6) x 10-^ 
in remarkable agreement with the CKM unitarity of the SM. This is illustrated in Fig. [TJ Next, we 
discuss the implications on BSM physics in a model-independent framework. 

3.1.1 Model independent constraints 

Each element Vij receives a universal, or channel independent, correction due to possible new physics 



corrections in muon decay, parameterized by in Eq. (2.8). Additionally, Vij receives channel- 
dependent BSM contributions, which are linear combinations of the e's defined in Eq. (2.7). Given 
the hierarchy |\4dp ^ iKsP, we discuss in detail only the BSM contributions to Vud, beginning with 
the O"*" — 7- 0^ nuclear transitions. 

From each 0"^ — )■ O"*" transition, working in the impulse approximation, one extracts the quantity 



|^tici|o+^.0+ ~ I Kid I ^ 



l + 2Re{eL + eR-e^) + c^+{Z) gsRees , (3.25) 
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Figure 1: Graphical representation of the current status of \ Vud\, \Vus\, and the corresponding 
CKM unitarity test. The horizontal band represents the constraint from decays, the thin 
vertical band the constraint from 0^ — ?■ 0^ nuclear decays, the oblique band the constraint 
from Kfj_2/'^ii2i whereas the ellipse is the la fit region. Figure reproduced with permission 
from Ref. pl^ . 



where the first correction reflects the BSM shift in the vector operator minus the shift in the Fermi 
constant extracted in muon decay. The second correction, proportional to Rees, arises from a non- 
vanishing Fierz interference term; this distorts the electron spectrum and therefore the phase space 
integals. The correction depends on the individual nuclear transitions, through [171 

4+{Z) = -2Vl - a-'Z^ hiQEclme) ^ I,(xo) = H x''" (xo - dx , (3.26) 

Io{QEc/me) Ji 

where Qec = Mi — M2 with rrie the electron mass. The transition strengths, or ft values, after the 
application of transition-dependent radiative corrections and isospin-symmetry-breaking corrections to 
the nuclear matrix elements, are remarkably constant with the Z of the daughter nucleus, supporting 
the Conserved- Vector-Current (CVC) "hypothesis", though CVC is simply a consequence of the SM. 
Moreover, CVC is also tested by studies of nuclear mirror transitions, albeit with lesser precision |1U9] . 
The computation of the isospin-symmetry-breaking corrections, noting e.g. Ref. [85], has been criti- 
cized [891 ^0], but the employed procedure has been experimentally validated in a system for which 
the correction is particularly large [120] . For a discussion of the various methods to compute isosp in- 
breaking effects, see Ref. ^1] . From the constancy of the corrected ft values with the Z of the daughter 
nucleus. Hardy and Towner [SSj obtain the combined bound 

- 1.0 X 10"^ < gs Re < 3.2 x 10"^ (90% C.L.) . (3.27) 

It is the most stringent bound on scalar interactions from low-energy probes. Moreover, the Ackm 
constraint implies 

Re(eL + eR-e^) <5x 10"^ (90% C.L.) . (3.28) 
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This is one of the strongest precision constraints on new physics. It corresponds to effective scales 
A > 11 TeV [IS]. The resulting constraints on the weak-scale couplings otj [ISl HH ETj are comparable 
to that obtained from Z-pole experiments, and are stronger than the ones obtained from a{e^e~ — ?■ qq) 
at LEP. 

In principle, neutron decay allows for an extraction of V ud free of nuclear structure uncertainties. 
Assuming no BSM effects one has = [4908.7(1.9)s/(r„(l + 3c/^)]i/2 ^ ^ ^n extraction of Vud 
competitive with nuclear decays requires Sqa/qa ~ 0.025% and (5r„ ~ 0.35s {5Tn/Tn = 0.04%). In turn, 
a determination of qa at the required level necessitates a measurement of the beta asymmetry A at the 



0.1% level. The expression for \V. 



ud I 



n^peu 



reads |17] 



\Vud? 


= \Vud\ 




n^peu 



l+2Re(ei + eK-e^) + 



3A2 



gsRees-l2\gTReeT 



6A2 



/oK) 1 + 3A2' 



(3.29) 

where A = Qa/ Qv^ a^o = Eo/me, with £"0 the electron endpoint energy. The constant c is a certain 0{1) 
number that depends on the specific experimental analysis used to extract A from measurements of 
the beta asymmetry A, presuming the presence of spectrum contaminations due to scalar and tensor 
operators. A different correction would appear if one extracted A from a different observable, such as a 
(see the discussion in Ref. |i7]). 



3.2 Lepton Universality 

The ratio Rt, = T{n ^ eh'['y])/T{n — )■ yuz^[7]) is helicity-suppressed in the SM, due to the V — A structure 
of charged current couplings. It is therefore a sensitive probe of all SM extensions that induce axial and 
especially pseudoscalar currents, as well as of non-universal corrections to the charged current lepton 
couplings. The quantity R-^ can be predicted very precisely in the SM because the leading hadronic 
input, namely, the pion decay constant F^, cancels in the ratio. Once one includes electroweak radiative 
corrections, hadronic structure effects do appear, and the SM prediction can be organized within the 
ChPT power counting as follows: 



SM 



Ae2p2 + Ae2p4 + ... 



mf 



mz 



mf 



(3.30) 



(3.31) 



The leading electromagnetic correction Ae2p2 corresponds to the point-like approximation for the pion |121l 
[122] • The NNLO (two-loop) correction A^2p4 has been calculated within ChPT in Refs. [1231 HM]- The 
two-loop effective theory results have been complemented by a la.Tge-Nc calculation of an associated 
counterterm and by summation of leading logarithms a" In" (m^/mg) [122] giving |123[ 1124] 



R^ = (1.2352 ±0.0001) x 10" 



(3.32) 



The central value of i?^ is in agreement with the results of previous calculations |122l 1125] , pushing the 
theoretical uncertainty below the 0.1 per-mille level. 



3.2.1 Model independent constraints 

The ratio R^ = T{7i — >■ ei/[7])/r(7r — )■ /iz/[7]) probes more than the effective low-energy pseudoscalar 
coupling ep defined earlier as the coefficient of the operator e(l — 75)t'e • ^750?. In fact, since (i) Rj, 
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is defined as the ratio of electron-to-muon decay and (ii) tlie neutrino flavor is not observed in either 
decay, this observable is sensitive to the whole set of parameters e'if and e'^^ defined by 



eff 



D 



G p 

71 



ud 



a/3 



1 - 75)z//3 ■ M75ci + ep^ e„(l + 75)z//3 ■ w^^^d 



(3.33) 



where a G {e, /i} refers to the flavor of the charged lepton and [3 G {e, /x, r} refers to the neutrino 
flavor. One generically expects SM extensions to generate non-diagonal components in e^sT- 
new notation the previously defined pseudoscalar, scalar, and tensor couplings read ep^T = ^^pst- 
is important to note that only ep and e^p can interfere with the SM amplitudes, while the remaining 
and each enter as an absolute square in the numerator and denominator of R^^. In summary, 
allowing for non-standard axial and pseudoscalar interactions and factoring out the SM prediction for 
i?^, one can write: |17] 



R-K _ 




2 

+ 


Bp ^e/J- 
me P 


2 

+ 


Bo ,eT 


2 


Bo ~ea 
m,^P 


2" 






^ + ^L ^R ^P 


2 

+ 


Bp ^^e 
P 


2 

+ 


Bo_^f^r 
rrifi P 


2 

+ E. 


Bo -fJ-a 
m^^P 



1 + A 



e/fj. y 



(3.34) 



where we note that the BSM couplings can be complex. In the above equation the factors of Bo/rrie^nep 
represent the ratio of the new-physics amplitude over the SM amplitude. The latter is proportional 
to the charged-lepton mass due to angular-momentum conservation arguments, while the former is 
proportional to (0|'U75(i|7r), characterized by the scale- and scheme-dependent parametei]^ 



Ml 



(3.35) 



Since B^^{fi = 1 GeV) = 1.85 GeV (using the PDG |126] central values for the light quark masses) and 
consequently BQ/rrie = 3.62 x 10^, i?^ has enhanced sensitivity to ep^, and one needs to keep quadratic 
terms in these new physics coefficients. We discuss bounds on = — ep and separately. 
First, setting = 0, inspection of Eq. (3.34) reveals that if the new-physics couplings respect e|?/me = 
ep"/m^, then = 1, and there are no constraints on these couplings. On the other h and, i f 

the effective couplings are all of similar size, one can neglect the entire denominator in Eq. (3.34), 
as it is suppressed with respect to the numerator by powers of me/m^. We will assume the second 
scenario operates. In this case the constraint in Eq. (3.34) forces the couplings ep , 



, ep , e|? to live in 



a spherical shell of radius uiJEq^RT^ /Rf^ ^ 2.75 x 10"^ centered at Re(e|,") = mjBo ^ 2.75 x 10-^ 
= ep = e^p = 0. The thickness of the shell is numerically 1.38 x 10~^ and is determined by the 
current combined uncertainty in i?f p [1271 ESHl [129] and Rf^ ^M- RT^/Rf^ = 0.996(5) (90% 
C.L.). This is illustrated in Fig. [2| where we plot the allowed region in the two-dimensional plane given 
by Re(ep) and a generic "wrong-flavor" coupling denoted by ep {x ^ e) — or this can be Im(ep) or 
any of the real or imaginary parts of e^ and e|?. Note that the allowed region is given by the thickness 
of the curve in the figure, thus enforcing a strong correlation between ep and ep . Since e^ and others 
of that sort are essentially unconstrained by other measurements, though we expect they can be of 
C(10^^), we can neglect all of the couplings but one to obtain a bound on that one. The resulting 
bounds using R.,, at 90%-C.L. are 



- 1.4 X 10"^ < Re (e^^) < 5.5 x 10-^ or - 2.75 x 10"^ < Im (e^') < 2.75 x 10"^ 



(3.36) 



'Note that the scale and scheme dependence of Bq{^) is compensated in physical quantities by the scale and scheme 



dependence of the Wilson coefficients e 



P ■ 
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Figure 2: Illustration of the allowed region in the two-dimensional plane Re(ep)-Re(ep ) 
(with X e) determined by Rt^, which is given by an annulus of thickness 1.38 x 10~^. In 
the absence of information on Re(e^), the 90 % CL. bound on Re(ep) is —1.4 x 10^^ < 
Re(e^") < 5.5 x 10"^. Note that Im(e^"), Re(e^"), and Im(e^°) are subject to the same bound 
as Re(ep ). Figure adapted from Ref. [Hj. 



Note that Re(ep) and Im(ep ), as well as Re(e^) and Im(ef?), are all subject to the same bound as 
Im(ep). Our results are in qualitative agreement with the findings of Refs. j58| [27]. 

Alternatively, neglecting the pseudoscalar couplings, one obtains the following combined limit on 
the axial combination of new couplings: 



- 4.5 X lO"'' < Re (eX - e''/) < 0.5 x 10"' . (3.37) 

Finally, we discuss how R.,^ is also sensitive to non-standard scalar and tensor couplings. As originally 
discussed in Refs. [571 EHl |59] , the pseudoscalar coupling ep can be radiatively generated starting from 



nonzero es,T- Hence, the stringent constraint in Eq. (3.36) puts constraints on the same e^^^ that can 
be probed in beta decays. The physics of the effect is this: once the scalar, pseudoscalar, and tensor 
operators are generated by some non-standard physics at the matching scale A, electroweak radiative 
corrections induce mixing among these three operators. Thus even if ep(A) vanishes at the matching 
scale, known SM physics generates a nonzero tp{^) at some lower energy scale /x via loop diagrams. 
The general form of the constraint can be worked out by using the three-operator mixing results from 
Ref. [SH]. The leading-order result is 

(l + 7PP log -\ + ef (A) log - + (A) 7tp log - (3.38a) 

113 ttl 9 / , ^ 

— — ^ 1.3 X 10"^ (3.38b) 

72 TT 

:i 6.7 X 10"^ (3.38c) 

15 ai n 

- — — ^ -7.3 X 10"^ , (3.38d) 

2 TT 



ep^(/i) = 
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3 ^2 _^ 


7pp = 




4 TT 
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where ai = a/ cos^ 9w and 0^2 = sin^ 9w are the U(l) and SU(2) ^ weak couphngs, expressed in terms 
of the fine-structure constant and the weak mixing angle. Setting ep (A) = and neglecting the small 
0{a/n) fractional difference between es^ri-^) and the observable es,Tif^) at the low scale, the constraints 
on es and-eT using R^, at 90% C.L. read 

-1.4 X 10-^ ^ , , ^ , , 5.5 X 10"^ 

1 X < 75P Re €5 + 7TP Re eT < . . . , . , 3.39 

log(A//i) log(A//i) 

and 

2 75 X 10~^ 

|75P Im(e5) + 7TpIm(er)| < — . , (3.40) 

log(A//i) 

Assuming log(A//i) ~ 10, so that, e.g., A ~ lOTeV and ~ 1 GeV, and using the numerical values of 
lSP,TP, one finds that the individual constraints are at the level of |Re(es') | < 8 x 10~^, |Im(e5) | < 4x 10~^, 
|Re(eT)| < 10-^ and |Im(eT)| < 0.5 ■ 10^1 These bounds become logarithmically more stringent as the 
new-physics scale A grows. It is worth noting that analogous studies are also possible in kaon decays, 
and new results are expected from NA62 at CERN [25] and TREK at J-PARC [I30]. 

Constraints on is^x can be worked out similarly [SH], resulting in |Re(e5)| ^ 5 x 10~^, |Im(e5)| < 
2.5 X 10~^, |Re(eT)| < 0.6 x 10~^, and |Im(eT)| < 0.3 x 10~^, which together with ip are the strongest 
low-energy bounds on the e couphngs [37] . 



4 Decay correlations and non-{V — A) couplings 

Differential decay distributions in beta decays are very sensitive to the Lorentz structure of the under- 
lying weak interaction, thus enabling searches for small non-{V — A) components. Following Ref. [19], 
one writes the differential decay distribution in the nuclear decay P ^ D e~ v {e^v) as a function of 
electron (positron) energy and lepton directions as follows. 



Ef. E,, Ep E,, 



+ 



dEJVLj^l^ (27r)5^" ^ M E^ E,E^ V 

(4.41) 

where P and D represent the parent and daughter nuclei, (J)/ J represents the parent nucleus polariza- 
tion, if any, and Pe^u are the electron (positron) and antineutrino (neutrino) three-momenta. We have 
omitted the additional parity-conserving term which appears if J 7^ 1/2, as indicated by the ellipsis. 
The coefficient b is the Fierz interference term, a is the electron-antineutrino correlation, A is the beta 
asymmetry, B is the antineutrino asymmetry, and the coefficient D is T odd in that the associated 
triple product of vectors is motion-reversal odd. All these quantities contain combinations of the Lee- 
Yang effective coefficients as delineated in Ref. [19], as does ^, and are related to our parameters as per 



Eq. (2.17). Additional terms are present if one can measure the polarization of the emitted electron or 
positron [l9]. Note, too, that the various correlation coefficients become E^ dependent once corrections 
of radiative and recoil order are included. 

The decay correlations can be measured in neutron and nuclear decays, and substantial progress is 
expected in the next few years. In neutron decay, both cold and ultracold neutrons, implying distinct 
experimental techniques and hence entirely independent sources of systematic error, are used to measure 
these correlations. In the future we can expect experiments poised to take advantage of cold neutron 
beams of much greater intensity at the FRM-III (PERC) 0, the New Guide Hall at NIST [m], and the 
SNS (Nab) [11] . Concerning nuclear decays, the development of atomic trapping techniques has allowed 
the precise detection of daughter nucleus recoil momenta, which in turn permits bettered measurements 
of the electron-antineutrino correlation parameter a. 
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In the absence of radiative corrections, recoil corrections, and BSM contributions, the correlation 
coefficients a{Ee), A{Ee), and B{Eg) reduce to simple expressions, while h,D = vanish. For example, 
for a pure Gamow- Teller decay we have the prediction ogt = ~l/3, whereas for a pure Fermi transition 
we have ap = 1- For mixed Fermi-Gamow- Teller transitions there is also a precise prediction once the 
ratio of Fermi to Gamow- Teller strengths is known — and this can be determined from the lifetime. In 
the case of neutron decay, which is a mixed transition, one obtains: 



— > 



1- 

1 + 3A2 



A{Ee 



2A(1-A) 
1 + 3A2 



B{E, 



2A(1 + A) 
1 + 3A2 



(4.42) 



where A = QaIgv and the limiting value of B{Ee), e.g., is termed So- 
Going beyond the SM, the dependence of the correlations a, b, A, B, and D on the short- distance 
couplings ei and ii can be determined using their dependence on the couplings Ci ± C- given in Ref. |19] 
and the relations given in Eq. (2.17). The full expressions are quite complicated, but simplify con- 
siderably if one considers the leading linear corrections only. In regards to these, the salient points 
are: 



As mentioned previously, the right-handed coupling to linear order induces the shift A — ?■ A = 
A(l — 2ei{). In order to probe from correlation measurements, one needs to know A = Qa/ov 
independently; this can come from a LQCD calculation. 

The scalar and tensor couplings es^r appear at linear order only through the Fierz interference 
term b and the analogous term b^ in the antineutrino asymmetry parameter, where b^ is defined 
by B{Ee) = -Bo + byirte/Ee. Different nuclear transitions probe different combinations of the 
BSM couplings. For example, the Fierz term 6 in a pure Fermi or Gamow- Teller transition 
probes exclusively the scalar or tensor coupling, according to b-p = gsRe{es) and 6gt = 

VT 



±(87(7TRe(eT))/A, where j = y/l — a'^Z^ and the sign distinguishes emitters [SO]- Mixed 
transitions such as neutron decay probe a linear combination of scalar and tensor couplings. For 
neutron decay one has [50] : 
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1 + 3A2 
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1 + 3A2 



gsReies)X-AgTReieT)il + 2X) 



(4.43a) 
(4.43b) 



Measurements of the correlation coefficients a. A, and B always include contributions from the 
Fierz interference term b, and are therefore sensitive to es.T to linear order. This dependence 
arises because correlation measurements involve the construction of asymmetry ratios |132j . and 
the dependence on b does not cancel in the asymmetry denominators. For example, in order to 
isolate A[Ef;) one constructs the ratio 

, , , N:(E^-N(E^ 

where N±{Ee) are the spectra corresponding to events with J ■ > and J -pe <0, respectively, 
so that sensitivity to b does indeed appear through the denominator. In general, asymmetry 
measurements probe 

+ W + (AS. »....}. (4.45) 
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where the elhpsis denotes other possible corrections of radiative and recoil order whose appearence 
depend on the correlation considered. It is worth noting that simultaneous analysis of a{Eg) and 
y4(i?e), e.g., yields more powerful constraints on the underlying BSM contributions than either 
correlation alone |133] . 

The dependence on la couplings appears only to quadratic order, together with additional effects 
quadratic in the eps. In this section, we have considered beta-decay correlations only, though strong 
constraints on BSM couplings also come from the study of light meson decays. We now review con- 



straints on the various BSM couplings which appear in Eq. (2.17) 



4.1 Model-independent constraints on scalar and tensor couplings 

We now summarize the current best constraints on the scalar and tensor structures, and highlight 
prospects for future improvements. Currently, the most stringent constraint on the scalar coupling es 



arises from 0+ — )■ O"*" nuclear beta decays, as discussed in Section 3.1.1 , whereas the most stringent 
bound on the tensor effective coupling arises from the Dalitz-plot study of the radiative pion decay 
7r+ — ^ e^Vel- That is, an analysis of the Dalitz plot of this decay from the PIBETA collaboration |134j 
puts constraints on the product Re(eT/T) of the short-distance coupling er and the hadronic form factor 



fx defined by jl35] 

{-i{e,p)\ua^y-^r4\'n+) = - ^h {Pi^e^ - Puef,) , (4.46) 

where and are the photon four-momentum and polarization vector, respectively. The analysis 
of Ref. |135] . based on a large- A^c-inspired resonance-saturation model, provides fx = 0.24(4) at the 
renormalization scale = 1 GeV, with the parametric uncertainty induced by the uncertainty in the 
quark condensate. The 90%-C.L. experimental constraint —2.0 x 10"^ < /TRe(eT) < 2.6 x 10~^, 



when combined with the above estimate for fx evolved to 2 GeV implies 

- 1.1 X 10-=^ < Re(eT) < 1.36 x 10"^ (90% C.L.) . (4.47) 

This is the most stringent constraint on the tensor coupling from low-energy experiments. The next 
best constraints arise from measurements of nuclear beta decays [29j. 

Bounds on scalar and tensor interactions can be obtained from a number of observables in nuclear 
beta decays, other than 0^ — t- 0+ transitions. Although these bounds are currently not competitive, we 
summarize them here for completeness. The leading sensitivity to scalar and tensor operators appears 
through the Fierz interference term b. Significant constraints on b arise from the electron-polarization 
observables [49j as well as from measurements of A and a in Fermi, Gamow- Teller, and mixed transitions. 
Here is a summary of current bounds on es^T [29] : 

• The most stringent constraint from the beta asymmetry in pure Gamow- Teller transitions (Aqt) 
arises from ^°Co measurements and implies |136] 

-2.9x10"^ < grReier) < 1.5 x 10"^ (90% C.L.) . (4.48) 

Similar bounds can be obtained from measurements of Aqt in ^^^In decay [137]: —2.2 x 10^^ < 
^TRe(er) < 1.3 x 10"^ (90 % C.L.). 

• Measurements of the ratio Pf/-Pgt from the longitudinal polarization of the positron emitted in 
pure Fermi and Gamow- Teller transitions |138[ 1139] imply 

4 

-0.76x10"^ < gsRe{es) + -gTReier) < 1.0x10"^ (90% C.L.) . (4.49) 

A 

^"Note that there is a factor of 2 difference in the normahzation of the tensor couphng ex compared to what was used 
in Refs. [551ITM]. 
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• Preliminary results have been reported on the measurement of the longitudinal polarization of 
positrons emitted by polarized ^°^In nuclei |14Uj . The corresponding 90 % CL. sensitivity to 
tensor interactions, \gT Re(eT)| < 3.1 x 10~^, is quite promising although not yet competitive with 
radiative pion decay. 

• Finally, the beta-neutrino correlation a has been measured in a number of nuclear transitions |141t 
1142^ 1143^ 1144] . The resulting constraints on scalar and tensor interactions are summarized in Fig. 7 
of Ref. |141j . In terms of the coupling constants used here, the 90 % CL. combined bound on the 
tensor interaction reads \gT^Q{^T)\ < 5 x 10"'^, again not competitive with radiative pion decay. 

Future improvements can be expected from both neutron and nuclear decay measurements. In the 
case of neutrons [H [5], future measurements of the beta asymmetry A [121 [131 ttH US], the antineutrino 
asymmetry B |T3] , the electron-neutrino correlation a [T71 [THl HH] , and the Fierz interference term 
h [I?! 120] should exceed 10~^ precision. On the nuclear side, measurements of a and h in the pure 
Gamow- Teller decay of ^He [22] should also reach the 10~^ level in precision. In Fig. [3] we summarize 
the current constraints on and er (horizontal bands) and assess the impact of future measurements in 
both neutron and nuclear decays, where we assume \h,h^\ < 10~^ from neutron decay and |&gt| < 10~^ 
from nuclear decays at 90% CL. The future neutron constraints are represented by the diagonal bands, 
whereas the constraint from ^He is represented by the vertical ocher band in the plot. In Fig. [3] the left 
panel represents the constraints using quark model input for the scalar and tensor matrix elements gs,T, 
whereas the right panel uses our preferred LQCD estimates [17], which still have a 50% uncertainty in 
Qs- In the near future one can expect Sgs/gs to reach the 20% level from LQCD, thus increasing the 
constraining power of these measurements, so that the thickness of the bands will shrink. With current 
uncertainties on gs,T, measurements of 10~^-level precision will probe effective scales A5 > 7 TeV. 



4.2 Nuclear and neutron probes of recoil effects 

SM recoil effects can be tested to high accuracy via experiments involving correlations between final 
state electron/positrons and the polarization or alignment of the parent state, which have the form 

"^'^ 1 + MEe) P^^ + F{E,)A il-z ■ pj ■ - S^) (4.50) 



where P =< Jz> / J is the polarization and A=l — 3<J^>/ ( J( J + 1)) is the alignment. The beta 
correlation A receives both leading order and recoil corrections, whereas the alignment correlation F is 
purely a recoil order effect. However, both have been measured as functions of E^, and provide measures 
of the recoil form factors. In this regard, it is important to point out that the axial tensor form factor 
d{q^), which vanishes in the SM from isotopic spin invariance in the case of transitions between isotopic 
analog states such as in neutron beta decay or tritium decay, is i n g eneral nonvanishing. Indeed, it is 
in general comparable in size to the weak magnetism term b{q^Y^. In the case of mirror transitions 
the size of the axial tensor must be identical for electron and positron decays, and this is subject to 
test. A theoretical prediction for the weak magnetism term b in terms of the difference between parent 
and daughter magnetic moments exists for transitions between isotopic analogs, and in the case of 
mirror transitions the size of the weak magnetism term is given by the electromagnetic Ml width of 
the transition from the excited isotopic analog state of the daughter nucleus — another CVC test [146 



One expects corrections to this result which are linear in isospin-breaking, both in 6(0) itself and from 
possible second-class current contribution in the SM. In neutron beta decay, through study of a{Ee) and 
A{Ee), the second-class current contribution to d{0) can be determined independently of 6(0), implying 

^^Note that the vanishing of d(0) for transitions between isotopic analog states is a SM prediction and is violated by 
so-called second class currents which can arise if quarks have an additional quantum number ^145j. 
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Figure 3: Current and prospective 90% C.L. allowed regions in the Re(es')-Re(eT) plane 
implied by (i) the existing bounds on bp and vr — )■ ewy (horizontal (green) band); (ii) projected 
measurements of b and b^, in neutron decay (inner (red) and outer (blue) bow-tie shaped 
regions, respectively) at the 10~^ level; (iii) projected measurements of 6gt at the 10~^ level 
from ^He decays (vertical (ocher) band). Left panel: hadronic matrix elements taken in the 
ranges 0.25 < gs < 1-0, 0.6 < qt < 2.3 [27]. Right panel: scalar and tensor charges taken 
from LQCD, gs = 0.8(4) and qt = 1.05(35). Note that by reducing the uncertainty in gs the 
constraint on es' from bp becomes stronger, independent of any future neutron measurement. 
The effective couplings es,T are defined in the MS scheme at 2 GeV. Figure adapted from 
Ref. gZ]. 



that the CVC test can be made without assumptions regarding second-class currents |147j . Generally, 
since the form of the recoil effects to all observables has been given, together with electromagnetic 
corrections, these SM predictions are all testable [76] . 

An alternate route for testing the SM recoil predictions is to utilize decays where the daughter state 
is unstable and itself decays 

i) electromagnetically such as in the mirror transitions in the A = 20 system to the 2+ 1.63 MeV 
excited state of ^°Ne, which in turn decays via photon emission to the 0^ ground state [148] ; 

ii) strongly such as in the mirror transitions in the A = 8 system to the 2+ 2.90 MeV excited state 
of ^Be, which in turn decays via the emission of two alpha particles [149] . 

In the former case there exists a beta-gamma correlation 




(4.51) 



whereas in the latter there is a beta-alpha correlation 



^'"^ l + G{E,) (P^Y-S,)-^'^. (4.52) 



dE,dn,dn„ \ V E, / 3^2 / Mv* 



noting V* is the velocity of the alpha particle in the daughter rest frame. Here the decay correlation 
coefficient G is purely of recoil order, so that it is sensitive to recoil form factors. The form of G as well 
as of the radiative corrections have been calculated [76] . 
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4.3 Couplings involving right-handed neutrinos: ei^R^s.T.P 

Neglecting neutrino masses, all the Ip couplings enter the decay rates quadratically, i.e., as per oc \ep\'^. 
Detailed expressions of the contributions to neutron and nuclear beta decay correlation coefficients can 
be found in Ref. [50j, though one needs to re-express the Lee- Yang couplings in terms of the and 



using Eq. (2.17). The corresponding bounds can be obtained from the analysis of Ref. [29j, in particular 
from the fits to beta decay data allowing for non-zero e/,, only, implying |Re(eL -|- ej?)| < 6.4 x 10~^, 
|Re(eL - ej?|) < 5.8 x 10~^, and es, er only, implying \gsRe{es)\ < 5.5 x 10"^, |5fTRe(eT)| < 2.1 x 10"^ 
at 90% CL. Decay correlations set also bounds on the imaginary parts of these couplings. For example, 
using the measurements of a in the decay ^^Ar |143j and ^He |150i I151j . one gets |Im(e5')| < 0.17 and 
|Im(eT)| < 0.05 at 90% CL. We expect that these bounds can be improved by future more precise 
measurements. 



4.4 T-odd correlations 

Triple-product decay correlations can only be motion-reversal-odd and thus cannot be true tests of 
T-invariance |152] . As a consequence final-state interactions (FSI) can simulate a seeming T-odd corre- 
lation without a fundamental violation of T-invariance. In beta-decays the energy release is sufficiently 
small that the FSI are electromagnetic in nature, so that they are calculable with minimal hadronic 
ambiguity at accessible levels of precision |153[ I154[ I155j . Nevertheless, under an assumption of CPT 
invariance, such "T-odd" correlations are sensitive to new sources of CP violation, and thus to new 
physics, though an observation of such a correlation in excess of SM expectations would not allow one 
to conclude that T itself is violated |156] . 

T-odd correlations have been studied in kaon, neutron, hyperon |157[ I158j . and nuclear beta-decays. 
In decay, namely, — )■ fi~^i"y, the transverse muon polarization is studied, and the expected SM 
correlation is small [159^ 1160] with respect to that possible in models of new physics |161t 1159] . Existing 
experimental studies are consistent with no T- violation in this and related processes |162^ 1163] , but new 
results of greater sensitivity are expected from TREK at J- PARC |164l I165j . In ordinary beta-decay, 
a T-odd correlation is possible only if the initial state is polarized, or if the final-state polarization of 
one of the particles is observed. In the decay of polarized neutrons one can construct J ■ {pe x p,^), i.e., 
the D correlation, or J ■ [pe x ae), i.e., the R correlation, if the polarization of the emitted electron 
oc (Je is observed. Recently significant experimental efforts in regard to each of these correlations have 
been concluded. The emiT collaboration has presented the best limit on D in beta decay |166t 1167] , 
finding D„ = (—0.94 ± 1.89 ± 0.97) x 10~^ in neutron decay, a substantial improvement over earlier 
measurements [168^ 1169^ [25] . As for Dpsi, the 0{a) correction vanishes in the zero-recoil limit, and 
Dpsi ~ 10~^ |153[ 1170] . This calculation has been updated to employ the techniques of heavy-baryon 
chiral effective field theory by Ando et al. |155] : they reproduce the Callan-Treiman result in 0{aQ / Mi^) 
with Q ~ M„ — Mp — rrie and include the leading piece of the N^LO correction to realize Dpsi with an 
estimated accuracy of better than 1%. In terms of our non-standard couplings D in neutron decay can 
be written as 



Dbsm 



1 + 3A2 



4A Im(e/j) + Sgsgr Im (ese^ + e^e^) 



(4.53) 



Neglecting small quadratic effects in scalar and tensor couplings, the emiT limit translates to the 
90% CL. constraint —5 x 10~^ < Im(e/j) < 3 x 10"'*, where we have used A [HI] for A. For Im(e5,T) and 
lm{es,T) the analysis is more involved. D provides bounds on products such as Re(eT)Im(e5), etc., so 
that no bounds on the imaginary parts can be obtained unless non-zero real parts of the exotic couplings 
are assumed or discovered. 
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Recent results also exist for R. In the first neutron experiment both transverse components of the 
electron polarization are measured, to yield both R and N \171\ . The correlation N probes J ■ and is 
appreciably non-zero from FSI; the experimental measurement is consistent with the SM expectation, 
offering a sensitivity check of the setup |172] . They find Rn = 0.004 ± 0.012 ± 0.005, limiting the 
imaginary parts of scalar and tensor interactions via: 



Rv 



-8(7T(2A + l)Im(eT) - 2gs Xlm{es) . (4.54) 



l + 3An 

In comparison the measurement of R in ^Li decay yields R»u = 0.0009(22) |173] . which through 

Rsu = -l — lm{eT) (4.55) 

3 gA 

implies a limit —3.1 x 10~^ < Im(eT) < 1.8 x 10"'^ at 90% CL. after Rpsi has been removed. Conse- 
quently, the Rn limit is most usefully interpreted as a limit on Im(es), namely —0.15 < Im(e5) < 0.11 
at 90% C.L. Looser constraints (at the 20%-level) on the imaginary parts of the scalar and tensor 
couplings come from the measurement of a in 0^ — ?■ 0^ transitions |143t I142j and in ^He decay |144] , 
respectively — note Fig. 22 of Ref. |166] for a useful compilation. 

The possibility of constraining T-odd P-even couplings, including P decay parameters, via T-odd P- 
odd observables such as EDMs has been considered earlier in Ref. [174J • The basic idea is that through 
loop diagrams involving electroweak gauge boson exchanges, T-odd P-even interactions can generate 
T-odd P-odd interactions. The resulting bounds depend on the scale at which parity invariance is 
restored |175lll76j . In this context, the D correlation has been recently studied in Ref. |177] . considering 
how large D can be in light of constraints from electric dipole moment (EDM) searches. Focusing on 
the leading contribution to D (proportional to Im(eR)) the authors show that via a loop diagram the 
same phase contributes to the neutron and other EDMs. If this is the leading (or only) contribution 
to the neutron EDM, then one can conclude that the neutron EDM currently provides the strongest 
constraint on D, which is 10 — 10^ times stronger than current direct limits on D, depending on the 
model. Of course, the bounds on D can be weakened if other operators, other than the 4-fermion ones, 
contribute to the neutron EDM, and this is natural in many theories, interfering destructively with the 
contribution proportional to Im(e/{). The numerical evaluation of the effect of such operators becomes 
a model-dependent question, and one anticipates that the connection between dtoD can be weakened 
completely, though detailed investigation is warranted [1781 . 

In radiative /3 decay one can form a T-odd correlation from momenta alone, so that one probes new 
spin-independent sources of CP violation. A triple momentum correlation has been previously studied 
in — )■ TT^l^i'i'y [T79j, and its sensitivity to physics BSM considered |180] . In decay both elec- 
tromagnetic and strong, i.e., pion-mediated, radiative corrections can mimic the T-odd effect, but the 
electromagnetic-induced FSI are orders of magnitude larger |179l I18li I182j . Finally, a spin-independent 
T-odd correlation can be constructed in the radiative beta decay of neutron and nuclei |183j . propor- 
tional to ■ {pe X p^), offering the opportunity to study the imaginary part of the pseudo-Chern-Simons 
term |184j first found as a consequence of the baryon vector current anomaly and SU(2)lxU(1) gauge 
invariance at low energies |185[ I186[ I187j . In Ref. |183] . the effect of FSI on this new correlation have 
been computed, establishing the baseline for possible future searches of BSM CP-violating interactions. 



5 Collider limits on non-standard CC interactions 

The BSM interactions probed at low energy can also be directly probed at high-energy colliders. The 
collider signals, however, depend on whether the particles that generate the 4-fermion interactions are 
kinematically accessible at the collider energies. Model-independent statements can be made under 
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Figure 4: Projected joint 90% CL constraints on Re(e5') and Re(er) from future beta decays 
measurements and the LHC at ^/s = 14 TeV. The low-energy constraints correspond to 0.1% 
measurements of B, b in neutron decay and b in ^He decay, under two different scenarios for the 
lattice QCD uncertainties in gs,T- The LHC bounds are obtained by requiring less than 3 ev- 
produced signal events with: (i) m-r > 2.5 TeV and 10 fb~^ of integrated luminosity (solid, red 
ellipse); and (ii) uit > 4 TeV and 300 fb~^ (dashed, yellow ellipse). Cuts are chosen to reduce the 
expected leading background to be below 1 event. To obtain the projection it is assumed no events 
are found. Note that the effective couplings es,T are defined in the MS scheme at 2 GeV. Figure 
adapted from Ref. 



the assumption that the non-standard interactions remain point-like at TeV-scale energies; i.e., the 
mediators are not accessible at the LHC. In this case collider searches directly probe the various non- 
standard couplings ea and e^?, that contribute to the parton-level amplitude for pp — )• ei/ + X. Since the 
non-standard amplitudes do not interfere with the SM amplitude (except for terms proportional to ei), 
the LHC probes at the same level both the real and imaginary parts of the new couplings. In Ref. |17] 
bounds are derived on es^r by analyzing LHC data in the pp ^ eu + X channel at ^/s = 7 TeV and 
1 fb~^ integrated luminosity. In Ref. [37] the analysis has been extended to all non-standard charged- 
current couplings with integrated luminosity of 5 fb~^, using both the pp eu + X and pp — )■ e~^e~ +X 
channels through SU(2)i gauge invariance. 

The pp — )■ eu+X channel is directly related to beta decays, since the parton-level process is ud — )■ eu. 
In order to put bounds on the BSM couplings one uses the (cumulative) transverse mass distribu- 
tion |188] , noting that the transverse mass of the lepton pair is defined as = \J2E^E!^{1 — cos A^ei^- 
At high rriT the SM background falls off while the BSM interactions would produce events and thus 
increase their number. Similarly, for pp — )• e^e~ + X one uses the dilepton invariant mass distribu- 
tion |189] . dubbed mee, to constrain the presence of possible contact interactions. 

A comparison of the best bounds available for each interaction from low- and high-energy experi- 
ments is shown in Tables |2| [s] (for Re(eQ,) and Im(eQ,)) and|4| [s] (for Re(ea) and Im(ea)). Note that in 
these tables we report only direct bounds, leaving out bounds on the real and imaginary parts of es,T 
and es,T from Rj, as they can be evaded by cancellation. All of the tabulated results refer to a bound on 
the absolute value of the parameter unless a range is specified. The main points can be summarized as 
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follows [37] (see also Ref. |19Uj ). For the pseudo-scalar couplings ep and Ip the low-energy constraints 
from pion decay are at the 10~^ level, which are very hard to reach at the LHC in the near future. The 
same applies to the vector interactions eL,R, for which both the CKM-unitarity bound {Ke{eL,R)) and 
the emiT bound {\m{eii)) are at the 0.5 x 10-Mevel. 

For scalar and tensor interactions with left-handed neutrinos, low-energy experiments have tradi- 
tionally yielded stronger bounds on Re(e5) and Re(eT), but the current LHC bounds have caught up; 
and both probes are at the 10~^ and 10~^ level for Re(e5) and Re(eT), respectively. In the next few 
years we expect improvements in the bounds from both the LHC and low-energy experiments, through 
neutron [I5l [HI [161 HTj [121 EQ] and ^He decay p2J measurements at the 0.1% level and beyond. Pro- 
jected future bounds from both beta decays at the 10~^ level and the LHC on Re(e5_T) are shown in 
Figure |4j These results show that low-energy searches with 10~^ sensitivity would have unmatched 
constraining potential, even in the LHC era. Concerning the imaginary parts, for Im(eT) bounds from 
low-energy {R correlation in ^He) and the LHC are at the same level, while for Im(e5) the LHC bound 
is stronger than the one derived from the R correlation in neutron decay. 

Finally, for scalar and tensor interactions with right-handed neutrinos. Is and ey, the LHC bounds 
are also at the 10~^ and 10~^ level respectively, significantly better than the current and future low- 
energy limits for both Re(e5_T) and Im(es'^7-). To match the LHC bound one needs measurements of 
the electron-neutrino correlation a in Gamow- Teller transitions at the level of SacT/O'GT ~ 0.05 
Similar remarks apply to Ir, for which the LHC bound is 5 x 10~^, and no significant limit is available 
from low-energy probes. 

Table 2: Summary of 90% CL bounds (in units of 10"^) on the real parts of non-standard 
couplings obtained from low-energy and LHC searches (5 fb~^ at ^/s = 7 TeV). In order 
to deduce the low-energy bounds on the scalar and tensor couplings we used gs = 0.8(4) and 
gr = 1.05(35) |17]. Using gs = 1.08(28) |82] would lead to the stronger bound iRefe^)! < 
0.4 X 10~^. The couplings es,p,T are evaluated in the MS scheme at /i = 2 GeV. 
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6 Model constraints 

Our discussion has focused on model-independent constraints emergent from precise universality tests 
and correlation coefficient measurements, in terms of the effective couplings ea and e^. As we have 
mentioned, ultraviolet extensions of the SM will generate these non-standard effective couplings at 
some level, which are then functions of the model parameters. Therefore, the results we have presented 

^^Within a given NP model the ratio r(7r —> ev)/T{'K fiiy) is likely to produce the strongest bound not only on 
ep and ep, but also on £5.^ and es.T, through their loop-induced contribution. However, since these bounds are based 
on naturalness arguments, they could be circumvented by cancellations between different contributions, so that in a 
model-independent analysis the LHC offers the strongest constraint on €s,t- 
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Table 3: Summary of 90% CL bounds (in units of 10~^) on the imaginary parts of non- 
standard couplings ea obtained from low-energy (using liadronic input as per Table [2]) and 
LHC searches (5 fb~^ at i/i = 7 TeV). The couplings es,p,T are evaluated in the MS scheme 
at /i = 2 GeV. 
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Table 4: Summary of 90% CL bounds, in units of 10~^, on the real parts of the non-standard 
couplings ia obtained from low-energy and LHC searches (5 fb~^ at ^/s = 7 TeV). In order 
to deduce the low-energy bounds on the scalar and tensor couplings we used gs = 0.8(4) and 
gr = 1.05(35) |17]. Using gs = 1.08(28) ^ would lead to the stronger bound |Re(e5)| < 
6.9 X 10^^. The couplings es,p,T are evaluated in the MS scheme at /i = 2 GeV. 
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can be used to constrain the parameter space of any SM extension. Moreover, within a given SM 
extension, the low-energy effective couplings will show peculiar dependencies on the underlying model 
parameters, resulting in correlations among low-energy beta decay signatures and other observables. 
In this subsection we briefly illustrate these ideas by discussing non-standard CC interactions within 
the Left-Right Symmetric Model |19T| 1192] and the MSSM, for which detailed studies can be found in 
Refs. |193l 11941 11951 11961 1197] . For more extensive reviews of underlying models we refer to [27] . 

6.1 Left-Right Symmetric Model 

The Left-Right Symmetric Model |191] is based on an extended gauge group SU(2)iX SU(2)/jX U(l), in 
which in addition to the SM gauge assignments, the right-handed fermions transform as doublets under 
SU(2)ij. After spontaneous symmetry breaking, the charged gauge bosons Wl and Wr mix into light 
SM-like Wl, which is predominantly left-handed, and heavier W2, which is predominantly right-handed, 
states which mediate charged current processes: 

Wl = Wl cos C + W2 sin C Wr = -Wi sin C + VTa cos C, (6.56) 

with the mixing angle sm( ~ {v/vrY proportional to the ratio of the weak scale over the scale at 
which the SU(2)/j group is spontaneously broken, vr ~ Mw^, thus leading to the breaking of parity. To 
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Table 5: Summary of 90% CL bounds, in units of 10~^, on the imaginary parts of the non- 
standard couphngs la obtained from low-energy (using hadronic input as per Table [2]) and 
LHC searches (5 fb~^ at i/i = 7 TeV). The couplings es,p,T are evaluated in the MS scheme 
at /i = 2 GeV. 
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leading order in {v/vr)\ this model generates the following correlated CC non-standard couplings 

J 



ei = = , ei? = -C , cl = -C , e/? = , (6-57) 



with no dependence on the lepton mass — all other couplings are vanishing. It is clear then, that this 
minimal and manifestly left-right symmetric model predicts no appreciable deviations from the SM 
in lepton flavor universality tests. On the other hand, the model predicts Ackm = — 2C, so that the 
mixing angle, and therefore the scale of spontaneous parity breaking, is strongly constrained by Cabibbo 
universality tests. Stronger bounds still emerge from the K^-Ks mass difference [198] . though the role 
of long-distance contributions to the mass difference limit the severity of the constraint |199l 1200] . Given 
the strong bounds on ^, the model would also predict unobservably small effects in decay correlations 
sensitive to |eL,_Rp ~ C^- 

6.2 CC processes in the MSSM 

Within the MSSM with R-parity, the CC effective couplings ea are generated through loop diagrams 
(vertex corrections and box diagrams such as those depicted in Fig. [s]), resulting in expressions that are 



not nearly as simple as the ones in Eq. (6.57) 



The chirality flipping couplings esj- [201] require the presence of left-right mixing between sfermions 
running in the box diagrams of Fig. |5| which is proportional to the small Yukawa couplings or the 
trilinear soft "A" terms. Ref. |201j analyzed the phenomenological constraints on such mixing and 
determined the range of the allowed contributions to the weak decay coefficients h and 5, arguing 
that they may provide unique probes of left-right mixing in the first generation scalar fermion sector, 
provided a precision between 10~^ and 10~^ can be achieved. 

Concerning the universality tests, in Ref. |197j it was shown that the magnitude of the corrections 
^CKM and Ag/^ can be on the order of 10~^, which is consistent with precision electroweak tests and LHC 
direct searches for supersymmetric particles. The size of the universality violations is controlled by the 
splitting in the squark versus slepton spectra (Cabibbo universality) or in the selectron versus smuon 
spectra (lepton universality). Moreover, Ref. |197j showed that a comparison of the first row CKM 
unitarity tests with measurements of Re/fi. can provide unique probes of the spectrum of first generation 
squarks and first and second generation sleptons, as illustrated in Figure [6] and explained in the figure 
caption. As a consequence, universality tests will be powerful diagnostic tools if supersymmetric partners 
are discovered at the LHC. 

Finally, a discussion of the impact of Cabibbo and lepton universality tests within the R-parity 
violating MSSM can be found in Refs. [1951 M. 
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Figure 5: A subset of the Feynman diagrams giving rise to non-standard CC operators in the 
MSSM. 



7 Conclusions 



In this article we have reviewed the role of precision beta decays measurements in probing physics 
beyond the Standard Model in the LHC era. As for all precision tests, theoretical calculations of the 
SM amplitudes play a crucial role in setting the stage for BSM searches. Here we have tried to convey 
the flavor of the needed theoretical inputs, emphasizing the increasing role played by lattice QCD 
both in the meson sector, in regard to the determination of Vus, and in the nucleon sector, in probing 
non-standard scalar and tensor couplings through decay correlations. 

Concerning BSM physics, we have emphasized a model-independent EFT approach to beta decays, 
assuming that new physics is emergent at high energies, based on a quark-lepton level effective La- 
grangian. This approach has a dual benefit. On one hand, it allows an unambiguous comparison of 
the physics reach of probes involving different hadrons, such as pions and nucleons, and nuclei, limited 
only by our ability to compute the requisite matrix elements. Moreover, as stated at the start, in the 
absence of an emerging picture of new dynamics from collider searches, the EFT analysis is the first 
necessary step to establishing the motivation and significance of this set of low-energy probes. The 
current bounds on the real and imaginary parts of the non-standard couplings eL,R,s,p,T (involving left- 
handed neutrinos) and eL,R,s,p,T (involving right-handed neutrinos) are summarized in Tables |2| |3| |4| |5} 
The outlook is very positive: the effective couplings of all the BSM CC operators involving left-handed 
neutrinos are currently probed or will be soon probed in low-energy experiments at the level of 10"'^ or 
better. This corresponds to probing maximal BSM physics scales A ranging from 7 TeV (for scalar and 
tensor interactions) to 11 TeV (for vector interactions), to (9(100) TeV (for pseudoscalar interactions, 
assuming no cancellations and no mass or Yukawa suppressions). 

In all cases, the effective scale probed overlaps with the LHC reach: therefore, if new particles are 
found at the LHC, beta decays will play an important role in the "LHC inverse problem", i.e. in es- 
tablishing the properties of the new BSM dynamics. This is also explicitly illustrated in the case of the 
MSSM (see Section [6] and Fig. |6]). Moreover, if new BSM dynamics is above the LHC reach, sometimes 
termed the "nightmare scenario," one can analyze LHC data on the process pp ^ eu + X in terms of 
the same EFT used at low energy, modulo the known QCD running of the various couplings. Even 
in this pessimistic scenario, recent theoretical results [37] show that beta decay measurements at the 
10~^-level can be more sensitive than the LHC in probing non-standard CC interactions, noting once 
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Figure 6: Correlation between Ackm and Ag/^ in the MSSM. The red points (dark grey) arise 
from a generic parameter space scan. The green points (hght grey) arise after applying the con- 
straints from precision electroweak tests. The black points arise after applying the constraints 
from direct searches at the LHC. The three branches correspond to the following scenarios for the 
sfermion spectra: the vertical branch corresponds to light squarks, which are been largely ruled out 
by the LHC, and heavy sleptons; the right branch corresponds to light smuons and heavy selec- 
trons and squarks; the left branch corresponds to light selectrons and heavy smuons and squarks. 
Figure reprinted with permission from S. Bauman, J. Erler, M. Ramsey-Musolf, "Charged current 
universality and the MSSM", atXiv:1205.0035 [hep-ph] [197]. 



again Tables |2| [s} |4| [sj and Fig. |4| and that 10~^-level measurements would have unmatched sensitivity. 
These considerations illustrate the relevance of precision beta decay measurements throughout the LHC 
era. 
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